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Dust emission from young outflows: the case of L 11 57 
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Abstract. We present new high-sensitivity 1.3 mm bolometer observations of the young outflow L 1157. These data show 
that the continuum emission arises from four distinct components: a circumstellar disk, a protostellar envelope, an extended 
flattened envelope -the dense remnant of the molecular cloud in which the protostar was formed-, and the outflow itself, which 
represents ~20% of the total flux. The outflow emission exhibits two peaks that are coincident with the two strong shocks 
in the southern lobe of L 1157. We show that the mm continuum is dominated by thermal dust emission arising in the high 
velocity material. The spectral index derived from the new 1.3 mm data and 850 fj,m observations from Shirley et al. (2000), 
is ~5 in the outflow, significantly higher than in the protostellar envelope (~3.5). This can be explained by an important line 
contamination of the 850 fim map, and/or by different dust characteristics in the two regions, possibly smaller grains in the 
post-shocks regions of the outflow. Our observations show that bipolar outflows can present compact emission peaks which 
must not be misinterpreted as protostellar condensations when mapping star forming regions. 



1. Introduction 

Mm and submm mapping of the dust thermal continuum emis- 
sion has proven to be a powerful tool to study the density struc- 
ture of low-mass protostellar condensations at different evo- 
lutionary stages (e.g. Motte et al. 1998; Chandler & Richer 
2000; Hogerheijde & Sandell 2000; Shirley et al. 2000; Motte 
& Andre 2000). However, virtually all protostars are driving 
powerful molecular outflows (Bachiller 1996), which are much 
more extended than the protostellar envelope size (e.g. Andre 
et al. 1999). Strong interaction between the outflow and the en- 
velope is observed in a few cases (e.g. L 1527, Motte & Andre 
2000). But besides such local perturbations, there is increasing 
observational evidence that the outflow itself can be respon- 
sible for (sub)mm continuum emission also at long distances 
from the driving source (Shirley et al. 2000; Chini et al. 2001). 
Since such continuum emission can be a valuable tracer of the 
outflow material, its detailed origin and structure clearly de- 
serves further study. 

The outflow located in the L 1 157 molecular cloud, at a dis- 
tance of 440 pc, is one of the few which have been detected 
so far in the continuum submm and mm emission (Shirley et 
al. 2000; Chini et al. 2001). This outflow is driven by an ex- 
tremely young Class source (Gueth et al. 1997). The south- 
ern lobe consists of two large overlapping cavities, each one 
ended by a strong bow shock, which are misaligned due to the 
probable precession of the ejection axis (Gueth et al. 1996, 
1998; Bachiller et al. 2001). The two shocks exhibit a rich 
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molecular composition (Tafalla & Bachiller 1995; Zhang et 
al. 1995; Avery & Chiao 1996; Bachiller & Perez-Gutierrez 
1997) which makes L 1157 the prototype of chemically ac- 
tive outflows (Bachiller et al. 2001). In this letter, we report 
new high-sensitivity bolometer observations of L 1 157 at mm 
wavelengths, which, in combination with previously published 
maps, provides basic information on the continuum emitting 
regions. 

2. Observations 

The observations were carried out in 1997 February and 1999 
December, with the IRAM 30-m telescope. The detector was 
the MPIfR 1.3 mm bolometer array (MAMBO), which had 19 
channels in 1997 and 37 channels in 1999. Each channel has 
an equivalent bandwidth of ^^70 GHz. 8 individual on-the-fly 
maps were obtained, with a different scanning direction for 
each coverage, from nearly perpendicular to almost parallel 
to the flow axis. The antenna beam width was 12". Skydips 
were observed every one to two hours to monitor the atmo- 
spheric zenith optical depth, which was found to be 0. 1 to 0.4. 
Flux calibration was achieved by observing Uranus. The final 
flux accuracy is estimated to be ^30%. The data were reduced 
with the IRAM package (NIC), using the EKH restoration al- 
gorithm. Fig. la presents the resulting image. The rms at the 
map center is ^3.5 mJy/beam. 

Beltran et al. (2003) obtained IRAM Plateau de Bure in- 
terferometric observations of the 1.3 mm continuum emission 
of the L 1157 protostar. The bolometer map presented in this 
paper was used to derive the short-spacing information and 
thereby complement the interferometric dataset. The resulting 
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' al. (2003). First contour and contour increment are 10 mJy/beam 

00^0 (3o'). The angular resolution is 2". The dash circle indicates the 

Right Ascension (J2000) interferometer primary beam (22"). 



map, which covers only the central region of the bolometer im- 
age, is shown in Fig. lb. 

3. Results 

3. 1. Four continuum components 

The data presented in Fig. 1 show that the mm continuum emis- 
sion in LI 157 arises from four different components. Their 
properties are summarized in Table 1 . 

Central source. The interferometric observations (Fig. lb) 
reveal the existence of a compact (< 1") source, of ^^80 mJy, 
located at 20h39'"06!24, 68°02'15':6 (J2000). Obviously, this 
structure is closely associated to the protostar itself, and it is 
tempting to identify it with the thermal dust emission of a cir- 
cumstellar disk. We refer to Beltran et al. (2003) for a more 
detailed discussion. 

Envelope. The central source is surrounded by an extended 
emission, that originates from the infalling (Gueth et al. 1997; 
Mardones et al. 1997) protostellar envelope. On the bolometer 
map, the peak flux density is 500 mJy/beam, and the flux inte- 
grated in a circle of 40" diameter is 960 mJy. Assuming a typ- 
ical dust temperature of 30 K, the corresponding mass would 
be ~2.1 Mq (Table 1). At the resolution of these observations 
(12"), the envelope is smooth and symmetric (the contour at 
half intensity is nearly round and has a diameter of ~13", cf. 
Fig. la). However, the interferometric map (Fig. lb) shows that 
the perturbation of the envelope by the outflow affects the in- 
ternal structure of the envelope, which is no longer isotropic at 
that resolution, but extended along the flow axis. 

Extended envelope. The L1157 protosteUar envelope is 
surrounded by an even more diffuse emission extended roughly 
east-west, i.e. approximately perpendicular to the flow axis 
(Fig. la). As shown in Fig. 2, this feature is also seen in the 



Component 


Size 


Flux 


Temp. 


Mass 






(mJy) 


(K) 


(Mq) 


Central source 


<1" 


80 


30 


0.15 


Envelope 


40" 


960 


30 


2.1 


Extended envelope 


140" X 60" 


320 


13 


2.1 


Outflow 


40" X 100" 


335 


25-60 


0.3-0.9 



Table 1. Properties of the four components of the 1.3 mm 
continuum emission in L 1157. The observed integrated fluxes 
are converted into masses assuming optically thin thermal dust 
emission, a uniform temperature, and a mass opacity k = 
0.01 cm^g~^. The outflow emission has not been corrected for 
the molecular line contamination (see text). 



850 ^m continuum map obtained by Shirley et al. (2000), as 
well as in the lower resolution NH3 (1,1) map obtained by 
Bachiller et al. (1993). It is also detected in quiescent emis- 
sion of i^CO (Gueth et al. 1997), DCO+, and N2H+ (BachiUer 
et al. 2001). Very likely this structure is a dense remnant of the 
molecular cloud in which the protostar was formed. The size 
of this envelope in the mm map is about 140" x 60", its flux 
density is roughly constant, ~12mJy/beam, and its total flux 
is 320±50mJy (after removing the 960 mJy corresponding to 
the central source and its envelope). If the dust temperature is 
13 K (as suggested by the NH3 observations of Bachiller et al. 
1993), the mass of this extended structure is ~2.1 Mq. 

Outflow. Besides the emission associated with the protostar 
vicinity. Fig. la shows that there is also an emission extending 
toward the south-east. A comparison with the CO (2-1) map of 
Bachiller et al. (2001), which has the same angular resolution, 
clearly shows that this emission is associated to the outflow (see 
Fig. 2): the continuum emission presents two main structures 
(labeled Bl and B2), whose positions and shapes correspond 
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Fig. 2. Comparison of the 1.3 mm continuum emission map of L 1157 with maps of other gas and dust tracers, (a) NH3 (3,3) 
emission, fromTafalla & Bachiller (1995). (b) CO (2-1) emission, from Bachiller et al. (2001). (c) A 1.3 mm continuum emission, 
see Fig. 1. (d) A850 /xm continuum emission, from Shirley et al. (2000). (e) NH3 (1,1) emission, from Bachiller et al. (1993). 



very well to those of the brightest CO emission. This is fur- 
ther demonstrated by the morphological agreement that exists 
between the mm continuum map and the NH3 (3,3) emission 
(Fig. 2) and the high-resolution SiO (2-1) maps of Gueth et 
al. (1998). The NH3 and SiO emission trace shocks in molec- 
ular outflows. B 1 does also coincide with the HH 375 Herbig- 
Haro object (Chini et al. 2001). We thus conclude that signif- 
icant emission arises from the shocked regions of the L 1157 
outflow. The peak flux densities of Bl and B2 are ~50 and 
~25 mJy/beam, respectively, and their integrated fluxes are 245 
and 90 mJy. 

3.2. Origin of ttie outflow emission 

The broad bolometer bandwidth includes many molecular 
lines, which are quite strong in LI 157, and could thus cre- 
ate a spurious continuum emission. At the Bl position, the 
CO (2-1) line-integrated emission is 155 Kkm~^ (Bachiller et 
al. 2001), which corresponds to ~8 mJy/beam in the bolome- 
ter band, i.e. 16% of the detected emission. Other lines are 
weaker, but their sum can however be important: adding the in- 
tegrated intensities of all lines (but CO) detected by Bachiller 
& Perez-Gutierrez (1997) at Bl yields llOKkm-^ This puts 
the lines contribution to about 30% of the emission detected by 
the bolometer. This value could possibly be somewhat larger, 
since the Une survey was not complete and not covering the full 
bolometer band. 

Fig. 2 also shows the SCUBA 850 fim map obtained by 
Shirley et al. (2000). This map looks very similar to the 
1 .3 mm observations. From Hirano & Taniguchi (2001), we can 
roughly estimate the CO (3-2) Une-integrated emission to be 
200 K km~ ^ at the B 1 position, which corresponds to ~20% of 
the bolometer measurement. Other lines must contribute and, 
as in the 1.3 mm case, could match the CO intensity and thus 
raise the line contribution to about 40%. 

We conclude that the lines contributions to the 1 .3 mm and 
850 fim bolometer maps are significant, but are unlikely to ac- 
count for the whole emission, which is thus dominated by con- 
tinuum emission. A comparison of both maps reveals a spectral 
index in the range 3 to 5.5 (see below), which is typical for dust 
thermal emission, and thus suggests that any free-free emission 



is negligible in these images. Dust is actually thought to play 
a key role in the chemical processes taking place in shocked 
outflows (e.g. Schilke et al. 1997). 

It is difficult to estimate the temperature of the dust in 
the outflow: the IRAS measurements (that would be needed 
to build the SED and then derive a temperature) do not have 
the angular resolution required to resolve the different compo- 
nents contributing to the L 1 157 emission. Tafalla & Bachiller 
(1995) used multiline NH3 analysis to derive outflowing gas 
temperatures of 60 to 80 K. These values are upper limits to 
the temperature of the dust, which has a very short post-shock 
cooling time (e.g. Clark et al. 1986). Chini et al. (2001) sug- 
gested a dust temperature of 27 K for HH 375 (Bl). From the 
measured mm fluxes, assuming optically thin emission and a 
temperature between 25 and 60 K, we estimate that the mass of 
the blueshifted shocked regions is ~0.3 to 0.9 M©. 

3.3. Spectral index 

Fig. 3 presents a map of the spectral index 
a computed between 850 //m and 1.3 mm: 
Q: = log(S'850/im/'S'i.3mm)/log(1300/850). This image shows 
systematic variations across the mapped region. The lowest 
values of the index, in the range 3 to 4, are obtained around the 
protostellar condensation. They are in good agreement with 
the average value a = 3.4±0.3 found by Shirley et al. (2000) 
in a sample of Class sources. The highest values are in the 
range from 4.5 to 5.5, and are measured towards the outflow 
region. 

Possible artifacts. A calibration error in one or both maps 
would translate into an incorrect scale for the spectral indexes, 
but would not introduce a position-dependent effect. However, 
line contamination of the bolometer measurements in the out- 
flow region can affect the apparent spectral index by intro- 
ducing an offset to the actual index of the continuum emis- 
sion. Simple algebra shows that the values discussed in the 
previous section result in a being increased by 0.35 at Bl, 
the position where the line contamination reaches its highest 
value. The variations of a at other positions are expected to 
be smaller. However, uncertainties on the line contamination, 
as well as possible calibration errors, make that this estimate of 
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the a variation is only a rough approximation. In fact, the offset 
could be smaller (or even negative), if e.g. line contamination 
at 1.3 mm was underestimated. Inversely, if the line contribu- 
tion at 850 iJ.m is higher, the offset can be larger: in order to 
introduce the observed 1.5 difference between the outflow and 
the protostellar envelope, a line/continuum ratio at 850 /im of 
~2 is required. This would imply that the 850 fim map is es- 
sentiaUy Une emission, a possibility which -from the fluxes of 
the lines observed in the band- seems unlikely. 

Dust properties. Besides the possible effects of the line 
contamination, a significant fraction of the variation in the 
spectral index has to be traced back to intrinsic properties of the 
continuum emission. For optically thin thermal dust emission, 
in the Rayleigh- Jeans limit, a = 2 + /3, where /3 is the expo- 
nent in the dust opacity law k oc u^. Hence, our observations 
would suggest that /? varies substantially from ~1.5 around 
the protostar to ~3 in the shocked regions^ The frequency 
variation of dust opacity is not well known, and depends on 
dust characteristics (see e.g. Visser et al. 1998, and references 
therein). Interestingly, a larger /3 - as suggested in the flow - 
points toward smaller dust grains (e.g. Kriigel & Siebenmorgen 
1994) and/or different chemical properties (e.g. Pollack et al. 
1994). Since the mm continuum emission in L 11 57 arises from 
shocked regions, as shown by the morphological coincidences 
with shock tracers such as NH3 and SiO, this result would be 
consistent with models predicting mantle evaporation (hence 
the rich chemistry observed in the shocks of L 1 157) and even 
grain destruction (which may be the only way to release siUcon 
and produce strong SiO emission, Schilke et al. 1997). 



4. Conclusion 

The observations presented in this letter demonstrate that the 
mm continuum emission in L1157 originates from four dis- 
tinct components: the circumstellar disk, the infalhng enve- 
lope, the dense remnant of the cloud in which the protostar has 
formed, and the outflow. The emission from the outflow repre- 
sents ~20% of the total flux at 1.3 mm, and is dominated by 
dust thermal emission. 

The spectral index computed between 1 .3 mm and 850 /um 
shows a significant difference between the protostellar enve- 
lope and the outflow emission. This could be related to an im- 
portant Une contamination of the 850 /xm map, and/or to varia- 
tions of the dust opacity law index. In the latter case, this may 
be an indication of different dust properties, possibly of smaller 
grains in the outflow. Observations of additional bipolar flows 
are needed to check if such variations in the spectral index are 
a common phenomenon, and to investigate their origin. 



' Note however that the factor hi//k is ~11 K at 1.3 mm and ~17 K 
at 850 nm, so the Rayleigh- Jeans approximation {T^lwlk) fails 
e.g. in the extended envelope, where Tdust~13K (if the dust is at 
the temperature indicated by the NH3 observations). In that case, 
Q = 2 + /3 + 7(rdust), where 7 can easily be derived from the Planck 
function (see e.g. Visser et al. 1998). At 13 K, 7 ~ —0.6, but this cor- 
rection is only 0.2 at 30 K, and < 0.1 for Tdust > 60 K. Hence, this 
temperature effect cannot explain the observed variations of a across 
the L 1 157 area, and one has to invoke variations of 
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Fig. 3. Grey scale map of the spectral index around the 
blueshifted lobe of LI 157, computed between 850 /itm and 
1.3 mm. The value of the spectral index is indicated by the scale 
shown on the right. Three contours of the 1 .3 mm emission map 
(Fig. 1) are also drawn for orientation. Before computing the 
index, we checked that both images agreed in position within 
~2", and smoothed them to the same resolution (20"). 

Finally, we caution that large-scale mm and submm contin- 
uum mapping of molecular clouds (e.g. Motte et al. 1998) or of 
the vicinity of protostellar objects (e.g. Shirley et al. 2000) can 
potentially be confused by the presence of outflow emission, 
such as that from LI 157. Some emission knots could easily 
be misinterpreted as, e.g., starless clumps. The observations re- 
ported here demonstrate that emission from bipolar outflows 
can be important, and must be taken into account when inter- 
preting the structure of mm maps of star forming regions. 
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